On the origin of the electric carrier concentration in graphite 
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We investigate the dependence of the electrical resistivity and magnetoresistance of single crys- 
talline micrometer-sized graphite samples of a few tens of nanometers thick on the defect concen- 
tration produced by irradiation at low fluences. We show that the carrier density of graphite n 
is extremely sensitive to the induced defects for concentrations as low as ~ 0.1 ppm and follows 
n ~ l/Ry with Rv the distance between defects in the graphene plane. These and Shubnikov-de 
Haas oscillations results indicate that at least a relevant part of the carrier densities measured in 
graphite is not intrinsic. 

PACS numbers: 73.90.+f,61.80.-x,81.05.Uw 



The electronic properties of ideal graphite are actually 
not well known simply because defect-free graphite sam- 
ples do not exist. In the last fifty years scientists flooded 
the literature with reports on different kinds of electronic: 
measurements on graphite samples, providing evidence 
for carrier (electron plus hole) densities per graphene 
layer at low temperatures no ~ 10 10 . . . 10 12 cm -2 , see 
e.g. Refs. [H, 0, 0. Taking into account that: (1) an 
exhaustive experience accumulated in gapless semicon- 
ductors - whose density of states should be similar to 
its counterpart in a semimetal - indicates so far that the 
measured no is most probably due to impurities Q , and 
(2) the expected sensitivity of no in graphite on lattice 
defects and impurities or adatoms @,la], a fundamental 
question remains unanswered, namely, how large is the 
intrinsic uq of ideal graphite? Why no is so important? 
Let us recapitulate some fundamental band structure 
theoretical results for the graphite structure Q. Two- 
dimensional (2D) calculations assuming a coupling 70 
between nearest in-plane neighbors C-atoms give a car- 
rier density (per C-atom) n(T) = (0.3 . . . 0.4)(fc s T/ 7o ) 2 
(70 ~ 3 eV and T is the temperature). Introducing a 
coupling (71 ~ 0.3 eV) between C-atoms of type a in ad- 
jacent planes one obtains n(T) = a(7i/7 2 )T + f7(T/ 7 o) 2 + 
c(T 3 /7 2 7i) + . . . (a, 6, c, . . . are numerical constants). In 
both cases n(T — > 0) — > 0. Neither in single layer 
graphene nor in graphite such T— dependences were ever 
reported Q, i.e. a large density background uq was al- 
ways measured and assumed as "intrinsic" without tak- 
ing care of any influence from lattice defects (including 
edge effects [7]) or impurities. To fit experimental data 
and obtain a finite Fermi energy Ep, up to seven free 
parameters were introduced in the past, whereas in the 
simplest case Ep oc 72 0,0- 

Clearly, any evidence that speaks against an intrinsic 
origin of - even a part of - the measured no in graphite 
samples would cast doubts on the relevance of related 
electronic band structure parameters obtained in the past 
and will help significantly to clarify observed transport 
phenomena. As in the case of gapless semiconductors |4j 
this requires a formidable experimental task. For exam- 
ple, to prove that the measured uq = 2 x 10 8 cm -2 in 



Ref. [9| is due to vacancies/interstitials requires a vacancy 
resolution better than 0.05 ppm. Although nowadays 
the concentration of impurities in graphite can be mea- 
sured with ~ 0.1 ppm resolution there is no experimental 
method that allows us to determine with such a precision 
the number of vacancies or C-interstitials. In spite of that 
and because of this situation we would like to start the 
discussion on the origin of n$ postulating that at least 
part of it cannot be intrinsic. The studies presented here 
provide answers to: (1) Can a single vacancy/interstitial 
provide ~one carrier into the conduction band even if 
they are several hundreds of nm apart (ppm concentra- 
tion)? This is a relevant issue specially because we ex- 
pect that the Fermi wavelength in graphite Xf > 1 M m 
[id ] . (2) Can the resistivity of graphite change with such 
small defect concentrations? (3) How reliable are band 
structure parameters of graphite obtained from the field- 
induced quantum oscillations in the resistivity (or magne- 
tization)? (4) Why there is an apparent maximum value 
for no ~ 10 12 cm~ 2 in graphite samples? 

In this study we measured the change of the elec- 
trical resistance of thin crystalline graphite samples as 
a function of defect concentrations between ~ 0.1 to 
~ 10 3 ppm. To do this we irradiated three ~ 60 nm thick 
and tens square micrometer samples under ambient con- 
ditions with a focused proton microbeam of 2.25 MeV en- 
ergy scanned over the samples. A photo of sample 1 can 
be seen in Fig. 1. Particle induced x-ray emission mea- 
surements were done in situ and revealed a total concen- 
tration of ~ 20 /xg/g of non-magnetic imp urities except 
hydrogen of concentration 0.5 ± 0.3% [llj . The magne- 
toresistance of a fourth sample of size 11 x 2 x 0.015 /an 3 
was measured at three different parts, each of length 
~ 1.6 fim and irradiated with 30 keV Ga + ions. In 
this sample Shubnikov-de Haas (SdH) oscillations were 
measured at 4K before and after irradiation and also 
for the corresponding bulk sample. Further details on 
the electron-beam spectroscopy and Raman techniques 
and the exfoliation and ultrasonic procedures use to char- 
acterize and prepare the crystalline thin graphite flakes 
from bulk samples will be published elsewhere. 

Figure 2(a) shows the relative change of the resistance 
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FIG. 1: Temperature dependence of the resistance of sample 3 
in the as-prepared state and after proton irradiation with a 
fiuence of 9 x 10 13 cm -2 (continuous lines). Note the large 
change in R(T) after inducing only ~ 3 ppm vacancy density. 
The dashed lines are obtained assuming a Fermi energy Ef = 
Ea + ksT and a T~ 2 dependence for the mean free path, see 
text for details. The inset shows an optical microscope picture 
of sample 1 with gold electrodes on top. 



vs. time during and after irradiation of sample 3 at 
297 K. The curves are obtained at different initial relaxed 
states after application of a certain proton fiuence. When 
the beam starts to hit the sample we observe a clear de- 
crease in the resistance, whose amount depends on the 
fiuence used and on the initial sample state. Figure 2(b) 
shows the resistance change relative to the sample virgin 
state. In this figure t = s means the time at which the 
beam stops irradiating the sample. Remarkable is that 
for all samples in the virgin state we observe a decrease in 
the resistance of a few percent for induced defect density 
< 3 ppm (average defect distance Ry > 100 nm!) that 
remains after several hours after irradiation, i.e. in the 
relaxed state R(t > lh), see curves (1-3) in Fig. 2(b). In 
Fig. 3 we show this relative change for samples 2 and 3. 
A further increase of the fiuence increases the resistance 
in the relaxed state, see Fig. 3. The explanation for this 
behavior is that defects increase the carrier density n, as 
theoretically suggested @. Because defects also act as 
scattering centers, both the carrier mean free path l{Rv) 
and n(Ry) have to be taken into account. 

We assume graphite as a structure composed of 
weakly-coupled graphene sheets 



12J. Within a factor 



of two the initial value at 297 K before irradiation for 
the carrier density is ~ 6 x 10 10 cm -2 and for the 
mean free path li ~ 50 nm [9j. The smallness of I at 
297 K in comparison to the sample size allows us to use 
the Boltzmann-Drude semiclassical approach. This is im- 
portant because for our sample sizes and at T < 150 K 
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FIG. 2: (a) Relative change of the resistance R(t) — 
R(0)/R(0) measured vs. time during and after the proton 
beam hits sample 3. The different curves were obtained 
at different starting relaxed conditions after irradiation of 
3.1, 4.0, 4.3, 6.4, 7.2, 9.0 x 10 13 protons per cm -2 , correspond- 
ing to the curve numbers 7, 9, 10, 14, 16, 19, from top to bot- 
tom. The small oscillations in the resistance observed during 
irradiation are an artifact due to the overlapping of the pro- 
ton current and the ac current of the resistance bridge, (b) 
Change of resistance relative to its value in the virgin state 
Ro vs. time. The time scale is taken from the time at which 
the beam does not hit the sample anymore, i.e. the mini- 
mum in AR/R(0) in (a). Note the decrease of resistance with 
irradiation for the first three curves (1-3 from bottom) even 
in the relaxed states (after 1 hour). As in (a) the different 
curves are taken from the sample at different initial states ir- 
radiated with fluences (1,2, ...,19) = (1, 1.5, 2.0, 9) x 10 13 
protons/cm -2 . 



there is no straightforward theoretical approach that in- 
cludes ballistic and diffusive scattering that allows us to 
obtain in a simple way n(Ry) from the resistance. Hall 
effect measurements are not necessarily preferred to ob- 
tain n(Rv) since: - the Hall signal depends on at least six 
unknown parameters (n, I, m* for electrons and holes in- 
dependently) that change with defect concentration and 
T; - conventional multiband approaches appear to be in- 
adequate for graphite [9J; - added to these difficulties, the 
Hall signal of graphite can be anomalous at T < 150 K 

M. 



Following SRIM simulations [14| the produced defect 
concentration at a proton fiuence of 10 13 cm -2 would 
be ny ~ 10 9 cm~ 2 . Assuming that each defect in the 
graphene plane increases by one the carrier number, the 
increase in carrier density after irradiating such fiuence 
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FIG. 3: Relative change of resistance vs. fluence for sam- 
ple 2 (open symbols) and sample 3 (close symbols). The 
upper x-axis shows the corresponding scales as average de- 
fect distance within a graphene plane Rv and the defect den- 
sity in parts per million (ppm). The triangles represent the 
change of the resistance in the relaxed state after irradia- 
tion relative to the virgin state Ro, he. R(t > ih)/Ro — 1. 
The circles are obtained from Rmin/Ro — 1 and the squares 
Rrain/R(t > lh) — 1. The curves were obtained from Eq. (1) 
with the following parameters: ny = Q.1/R V ,U — 50 nm 
(long dashed); 0.1/7? 2 ,20 nm (dot); 1/7^,50 nm (contin- 
uous); 3/i?v,150 nm (red dash-dot); 3/Ry, 50 nm (dash- 
double dot). The short-dash curve was obtained assuming 
the usual 3D relation 7? oc 1/ln and with ray = 0.9/R v ,k = 
150 nm. Note that with this 3D relationship no minimum in 
the measured range is obtained within a broad variation of 
parameters. 
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or in terms of the related wave 
-l 



vector ky ~ (-nny) 1 / 2 = 5.6 x 10 4 cm -1 ~ O.lfci?- This 
simple estimate reveals that such small defect concentra- 
tions are relevant for the transport. 

At 297 K the produced defects by irradiation are 
metastable [15j, see Fig. 2. Therefore, we plot in Fig. 3 
the relative change of the resistance just at the end of the 
irradiation R rn in with respect to the virgin state Rq, i.e. 
Rmin/Ro — 1 (close and open circles in Fig. 3). Another 
possibility is to plot the relative change with respect to 
the resistance taken one hour after irradiation R re i, i.e. 
Rmin/Rrei — 1 (close and open squares in Fig. 3). Both 
ways minimize the influence of annealing effects and pro- 
vide a similar behavior. These relative changes indicate 
that the resistance reaches a minimum ~ 90% of its ini- 
tial value at Ry — 20 ... 30 nm, see Fig. 3. At higher 
fluences the resistance increases because the decrease of 
I starts to overwhelm the increase in n. 

A quantitative description of these data can be done 
taking into account the two dimensional resistivity [HI 
p = 2/e 2 v F N(Ep)Tp, where vp,Ep,Tp are the Fermi 
velocity, the Fermi energy and the scattering relax- 



FIG. 4: First derivative of the magnetoresistance measured 
at 4 K in two parts (a) and (b) of similar area 2 x 1.6 (im 2 
separated by 3 /im in sample 4 (0 — 20 fl/T range), (c) Sample 
part (a) after irradiation with a fluence of 5 x 10 11 Ga + -ions 
per cm 2 (20 — 80 Q/T range). The dashed curve was obtained 
for the bulk sample of size 2 x 1 x 0.2 mm 3 at 4 K. The 
magnetic field was applied normal to the graphene layers. 



ation time at Fermi energy. Using N(Ep) for clean 
graphene the expression for Ep = hvpkp and 
Tp = l/vp, one arrives at the simple expression p = 
(Tr/2) 1 / 2 (H/e 2 )l~ 1 n~ 1 / 2 . Furthermore, the carrier den- 
sity increases as n = rii + ny, where ny = l/Ry for one 
carrier per defect. Following Mathiessen's rule, the mean 
free path is given by l^ 1 = + ly , where U is the ini- 
tial value due to all scattering centers before irradiation 
and ly the mean free path due to the produced defects. 
The relative change of resistance can be written as 
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The solid curve shown in Fig. 3 is obtained with ly — 
1.15 x 10 6 [cm -1 ]i? v [cm 2 ] for U = 50 nm and ny = \/R v . 
Within logarithmic corrections, the obtained ly(Ry) 
function agrees quantitatively with that found in Rcf. 5. 
In Fig. 3 we show also other curves obtained using other 
values for and pre- factors for ny(Ry) as well as assum- 
ing the usual 3D relationship R oc 1/n instead of 1/n 1 / 2 . 
The comparison indicates that within a factor of two ny 
is indeed given by 1/R V (for Ry > 10 nm) and that the 
usual 3D relationship for R cannot describe the observed 
behavior within a reasonable range of parameters. 

The remarkable increase in R(T < 300 K) and the ob- 
served change in the temperature dependence of graphite 
after inducing only ~ 3 ppm defect density (Ry ~ 
100 nm) is mainly given by the decrease of Ep ~ 
Ep(0) + ksT oc yj n(T) with temperature. As shown in 
Ref. |9| Ep is basically determined by thermal electrons 
(note that E F ~ 330 K for n = 6 x 10 10 cm- 2 ) and its 
T-dependence overwhelms that of l(T). With E F (T), 
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h(T) oc T~ 2 Q and the parameters obtained from Fig. 3 
in Eq. (1) one can understand the observed temperature 
dependence of the resistance above ~ 100 K, see Fig. 1; at 
lower temperature the Boltzmann-Drude approach looses 
its validity. 

The band parameters of graphite were obtained mostly 
on macroscopic samples and usually from magneto- 
optical studies, SdH and de Haas- van Alphen oscillations, 
cyclotron resonance, etc. We doubt that in graphite sam- 
ples the defect density is negligible and therefore we ex- 
pect that the carrier density is neither small nor homo- 
geneously distributed. Within the 11 /Ltm length of sam- 
ple 4 we measured the magnetoresistance at 4 K and 
calculate its first derivative in different parts of similar 
area. The SdH oscillations depend on the sample posi- 
tion, see Fig. 4, indicating clearly inhomogeneities in the 
carrier concentration within micrometers in agreement 
with EFM results that revealed sub-micrometer domain- 
like carrier density distributions in graphite surfaces [13] • 

For the measured sample area that gives curve (a) in 
Fig. 4 and within experimental resolution there are no 
SdH oscillations up to a field B ~ 1.8 T in clear con- 
trast to the bulk sample, see Fig. 4. This fact can be 
understood assuming that in most of this sample part 
n o i$ 10 9 cm 2 . Then, the corresponding Fermi wave- 
length Xp > 0.8 «m is of the order of the sample size 
and larger than the cyclotron radius r c — m*vp/eB for 
B > 0.07 T assuming m* = 0.01m (m is the free elec- 
tron mass) . In this case we do not expect to observe any 
SdH oscillations. However, for B ~ 1.8 T and 2.8 T two 
maxima are observed. From the measured "period" P in 
1/B as well as from the first field at which the first maxi- 
mum appears we estimate the existence of domains of size 
< 2r c < 100 nm in which Xp < 50 nm, i.e. domains with 
n o <] 10 n cm" 2 within a matrix of much lower carrier 
concentration. This indicates that the description of the 
SdH oscillations in real graphite samples can be achieved 
only within the framework of inhomogeneous 2D systems 
lallSl, an issue rarely studied in the past. 

The selected Ga + irradiation produced an average de- 
fect concentration of ~ 10 12 cm" 2 , i.e. ~ 10 3 ppm (Ga 
implantation < 1 ppm) in the thin graphite sample, "ho- 
mogenizing" its carrier density distribution. After irra- 
diation the SdH oscillations are clearly observed for B > 
0.7 T, see Fig. 4. Their period 0.16 T" 1 < P< 0.23 T" 1 
is within the range found in literature [2j, |2fJ and indi- 
cates no ~ 3 x 10 11 cm -2 . 

Finally, if a relevant part of the reported carrier con- 
centration in graphite is due to defects, why does it ap- 
pear to saturate at no ~ 10 12 cm" 2 (~ 3 x 10" 4 carrier 
per C-atom) ? We note that such saturation is also ob- 
served in gapless semiconductors with increasing donor 
concentration 0]. We might therefore expect it when 
the average distance between defects is of the order of 
the range of modification of the electronic structure pro- 
duced by, e.g. a single vacancy, found experimentally to 



be ~ 3 nm [2l[ implying uq < 10 13 cm 2 . 

Concluding, the obtained results indicate that a con- 
centration of defects (or impurities) of ~ 0.2 ppm can 
generate a carrier density ~ 10 9 cm -2 affecting the trans- 
port properties. This is an extraordinary sensitivity tak- 
ing into account the large distances between defects these 
concentrations imply. Taking into account that, in best 
case, we have an impurity concentration < 20 ppm, ex- 
cept for hydrogen < 1%, plus an unknown concentration 
of vacancies and interstitials, we should doubt about the 
assumed "intrinsic" origin of the measured carrier con- 
centrations in graphite. The behavior of SdH oscillations 
in micrometer-sized graphite regions and their changes 
after introducing defects support the above statement 
and indicate that real graphite is composed by an in- 
homogeneous distribution of carrier density. 
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